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The fermionic Gaussian operator basis provides a representation for treating strongly cor- 
related fermion systems, as well as playing an important role in random matrix theory. 
We prove that a resolution of unity exists for any even distribution of eigenvalues over 
hermitian fermionic Gaussian operators in the nonstandard symmetry classes. This has 
some important consequences. It demonstrates a useful technique for constructing funda- 
mental mathematical identities in an exponentially complex Hilbert space. It also shows 
that, to obtain nontrivial results for random matrix canonical ensembles in the nonstandard 
symmetry classes, it is necessary to consider ensembles that are not even functions of the 
eigenvalues. We show that the same restriction does not apply to the standard Wigner- 
Dyson symmetry classes of random matrices. 
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I. INTRODUCTION 



The linear transformations of fermion operators that preserve anti-commutation relations is an 
important fundamental symmetry group. It is known to have a representation as an exponential of 
a general quadratic form in the fermion operators. By analogy with the Gaussian function in statis- 
tics, we term this a fermionic Gaussian operator. The group properties and integration measures of 
these operators have been extensively investigated in the mathematical physics literature They 
have been used to generate the generalized fermionic coherent states^"-. This group is also known 
to correspond to a general symmetry class in random matrix theorj^^i^^i^. 

The Gaussian fermion operators have other fundamental applications. They are an operator 
basis for all fermionic density operators. Both the well-known thermal states and the BCS states 
are examples of this. The Gaussian fermionic phase-space representation^^ril makes use of this 
property to generate a positive distribution over such operators. This is a complete phase-space 
representation for either fermionic and bosonic many body density operators. It has been used, 
for example, to computationally evaluate the ground state of the fermionic Hubbard model and 
to obtain an expression for the linear entropy of a density operator—. 

In this communication, we show that an expansion using hermitian fermionic Gaussian oper- 
ators provides a continuous positive-definite resolution of the fermionic identity operator. Our 
result holds for any integrable even distribution of eigenvalues. This resolution is complete, in that 
it generates an identity operator for the entire Hilbert space, not a subset. The proof of the resolu- 
tion of unity is given in terms of the polar coordinates of skew symmetric matrices. This shows that 
the fermionic Gaussian operators can be thought of as coherent operators. In terms of continuity 
and completeness, they have an analogous role to the bosonic coherent state projectors'*^. 

To understand this comparison, we recall that one of the most useful properties of the bosonic 
coherent states is their resolution of unity'^'^U. This resolution of the bosonic identity operator 
uses positive-definite coherent state projection operators, with a positive measure over complex 
amplitudes. It is the basis for such well-known methods as the Husimi Q-function^. Our res- 
olution of unity for the fermionic Gaussian operators has similar properties, with a normalized 
Gaussian operator replacing the coherent- state projector in the expansion of the unity. Just as in 
the boson case, a resolution of unity is the foundation for calculating other fermionic identities and 
phase-space representations. These will be treated in subsequent publications. 

An alternative approach by Gilmore^2j2^ and Perelomov^li^ defined coherent state representa- 
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tions for fermions using generalized fermionic coherent states. These are obtained as the action 
of a representation of a Lie group on an extremal state. This definition is related to the dynam- 
ical structure of the corresponding creation and annihilation operators, using U (N) Lie group 
methods2r2ii^a2^~i2. In this approach, an integral over the coherent state projectors appears, but it 
does not generate an identity operator for all fermion states, due to number conservation. Because 
the fermionic coherent states are defined only for a subset of the Hilbert space, it is not possible to 
obtain a continuous representation of unity. 

Phase space representations for fermions^ ^ can also be defined in terms of Grassmann coherent 
states^i^"-^, using anti-commuting Grassmann variables^^. As the Grassmann coherent states 
form a complete set, it is possible in this case to obtain a resolution of unity. However, the use of 
Grassmann variables means that the phase-space is defined over non-commuting variables. These 
have an exponentially large matrix representation in terms of standard complex variables. Conse- 
quently, they have no efficient direct computational representation without requiring exponentially 
large computational resources. In addition, while Grassmann coherent states are mathematically 
well-defined, they do not correspond to any physical state. 

There is a close relationship between our work and the theory of random matrices. The trans- 
formations we use correspond to the nonstandard symmetry classes of random matrices. In this 
context, an important consequence of our result is that it provides an exact solution to a random 
matrix distribution of fermionic canonical ensembles. We show that in the nonstandard symmetry 
classes, an even distribution over matrix eigenvalues reduces a general canonical ensemble to the 
identity operator. In other words, canonical ensembles in these symmetry classes are nontrivial 
only if the random matrix distributions are not even functions of the eigenvalue. Nonstandard 
classes of transformations can still have useful properties if the eigenvalue distributions are not 
even functions - for example, if there is a lower cutoff to the eigenvalue distribution. Similarly, 
our result clearly does not exclude the use of ensembles that are not canonical. 

This paper is organized as follows: In section HI] we define the general Gaussian operators. 
In section Hill we define the symmetry classes that we use. In section |IV] we give the resolu- 
tion of unity for the general Gaussian operators, and in section |V] we consider the case of the 
number-conserving Gaussian operators. Section |Vl] gives a summary of our results, conclusions 
and outlook. 

Additional results and integration identities are given in the Appendix. 
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II. FERMIONIC GAUSSIAN OPERATORS 



We start by reviewing the known properties of fermionic Gaussian operators. This will help to 
establish the notation, as well as to give results that will be used in later parts of the paper. We 
consider a fermionic system that can be decomposed into M spatial or internal modes. In some 
treatments of symmetry properties the space and spin indices are treated separately, but this is not 
essential to our results. We define a as a vector of M annihilation operators and as vector of M 
creation operators. As they are Fermi operators, di and aj obey the anticommutation relations: 

|a,-,a]| = dij, {di,dj} =0. (2.1) 

We define an extended vector of all 2M operators y = (a,a^), with an adjoint vector defined 
as y' = {a' ,a) — . . . , a|^, ai , . . . , j • In the remainder of this section, we will summarize 
properties of the fermionic Gaussian operators to be used later. 



A. Group properties 



Following the original work of Balian and Brezin^ , we consider general linear transformations 
of fermions. These are obtained from the fermionic Gaussian operator defined as an exponential 
of a general quadratic form in Fermi operators: 



Gb (R) = exp 



(2.2) 



where R is a 2M x 2M antisymmetric complex matrix. The group composition law for the general 
Gaussian operator is 

GB(R) = GB(R(i))Gfi(R(2)), (2.3) 
The value of R under this composition law is obtained on defining: 



H = aR, 



(2.4) 



where we have introduced a symmetric matrix: 

Om Im 

With this definition, the matrices follow a composition law for the matrix parameters given by: 

exp (H) = exp (h^^^^ exp (^H^^^ 



(2.5) 



(2.6) 



This is clearly a Lie group. From the composition law, it has an inverse and an identity, as well 
as being associative and differentiable. It is equivalent to the 2M x 2M complex orthogonal Lie 
group, which follows^ from the anti-symmetry of R. 

B. Hermitian sub-group 

The general Gaussian operator for fermions can also be usefully defined in terms of the H 
matrix. Since this is the form most commonly used in later work, we will use it here. With this 
definition, the general un-normalized fermionic Gaussian operator is: 



Here H is an 2M x 2M matrix, whose definition in terms of an antisymmetric matrix R implies 
that it has the decomposition 



where h and A are M x M matrices. Hence the quadratic term in the exponent can be rewritten as: 



We note that, as well as being useful for linear transformations, the same class of fermionic 
Gaussian operators can also be used to expand a general fermionic density matrix in a posi- 
tive phase-space distribution—"—. In general, this distribution includes both hermitian and non- 
hermitian operators. 

In this form it is clear that if we wish to restrict the Gaussian operators to be hermitian, we 
must require that h = h^, and A = — A^. This is the form we will use here to establish the reso- 
lution of the fermionic identity operator. With this hermitian restriction, the quadratic form has a 
clear physical identification. It is simply the Bogoliubov-de Gennes Hamiltonian obtained from 
linearizing the Hamiltonian for a superconductor. From this perspective, we see that the Gaussian 
operator can have intrinsic coherence properties that correspond to either a superconductor with 
A 7^ 0, or to a normal fluid with A = 0. Because it carries phase information, the matrix A plays a 
similar role in fermion physics to the bosonic coherent state amplitude, which appears in various 
approaches to laser and superfluid theory. 




(2.7) 




(2.8) 




(2.9) 
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C. Number-conserving case 



Next we consider a subset of the general Gaussian operators with A = 0. These are the number- 
conserving fermionic Gaussian operators. In this case, H corresponds to the fermionic Hamil- 
tonian of a non-interacting Fermi gas in an arbitrary spin-dependent potential. In general, these 
Gaussian operators form a subgroup of the complex linear group, GL{M,C). Balian and Brezin^ 

found that these operators can be rewritten in terms of a single M x M matrix: 

1 



Gn (h) = exp 



a^ha--Tr{h) 



(2.10) 



Using the matrix identity exp(Tr(h)) = det(exp(h)), the number-conserving Gaussian opera- 
tors can be written as: 



G^(h) 



1 



exp 



/det(exp(h)) 

The group composition law for these operators is given by 



a^ha 



(2.11) 



G^(h) = Gn (h'^'>) Gm [y^^) , (2.12) 

Here the composition law for h is obtained on matrix multiplication of the variables \x = e^. The 
reason for this is that multiplying the u'-'^ matrices is equivalent to multiplying the exp ^H^-'^ 
matrices defined in Eq. (|2.6I) . Hence, we see that: 



(2.13) 



We can equivalently express this Gaussian operator in terms of normally-ordered parameters, 
using the following general mathematical identity for an M-mode fermionic operator—: 



Therefore we obtain: 



exp 



GAr(h) 



a^ha 



e 2 



: exp 



a 



e'^-l 



a 



Trfhl 



: exp 



a' 



gh-l 



a 



(2.14) 



(2.15) 



Hence, we notice that if we use normal ordering the simplest parameterization is through the 
parameter u = e^, which gives a matrix-multiplication group composition law for these operators. 
The quadratic term in the exponent can also be written in a Hamiltonian-like form, in terms of a 
thermal or number-conserving Hamiltonian Hyy: 



H, 



N 



1 



^ i a^ha — ah^ 



(2.16) 



We recognize that to obtain a hermitian Gaussian operator, it is necessary to restrict h to the class 
of hermitian matrices. 
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III. SYMMETRY CLASSES AND GROUPS 



The study of general symmetry groups for many-body systems originates in the work of 
Wigner— 1^ and Dyson^i"^, who studied the energy levels of complex many body systems such as 
nuclei. In particular, Dyson^i^ classified many-body systems depending on the symmetry proper- 
ties of the Hamiltonian. He classified Hamiltonians according to their time-reversal and rotational 
invariance properties. 

In random matrix theory, this classification corresponds to three random matrix models: the 
Gaussian unitary ensemble (GUE), the Gaussian orthogonal ensemble (GOE), and the Gaussian 
symplectic ensemble (GSE)^'-^^. The unitary ensemble applies to general systems without any 
invariance under time reversal, the orthogonal ensemble is for systems with time-reversal and 
rotational invariance, while the symplectic ensemble is for systems with time-reversal invariance 
but no rotational invariancd^i^. 

Dyson's classification is based on number-conserving Hamiltonians. If we consider the more 
general class of non-number-conserving Hamiltonians, there are four additional symmetry classes 
that can be identified depending on their time reversal symmetry and spin-rotation invariance^-. 
These four symmetry classes are Class D, Class C, Class Dili and Class CI, which are related to 
Cartan's classification of symmetric spaces. Class D corresponds to systems with neither time- 
reversal symmetry nor spin-rotation symmetry. Class Dili corresponds to systems with time- 
reversal symmetry but no spin-rotation invariance. Class C is for systems with spin-rotation invari- 
ance and no time-reversal invariance, while Class CI is for systems with spin-rotation invariance 
and time-reversal invariance. 

A. Gaussian operators in nonstandard symmetry classes 

The nonstandard symmetry classes are the general symmetry classes of matrices that involve 
non-number-conserving Hamiltonians and particle-hole symmetry. Here we focus mainly on the 
Class D symmetry defined by Altland and Zimbauer^, although our results are applicable to any of 
the nonstandard symmetry classes. Class D corresponds to the most general nonstandard symme- 
try class, as it applies to cases which do not have time-reversal or spin-rotation symmetry. In this 
case the only required properties of the matrices h and A, defined in Eq. (12.81) . are the hermiticity 
of h and skew-symmetry of A. The 2M x 2M H matrix is therefore hermitian, as discussed above. 
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We start by summarizing group theoretic results of Altland and Zirnbauer, then extend these to 
obtain new results for distributions over canonical ensembles. 

Hermitian matrices do not themselves form a group because they do not close under com- 
mutation. However, the anti-hermitian matrices form a Lie group, so it is useful to define an 
anti-hermitian matrix X, as X = /H. Hence the conditions on h and A can be expressed in terms 
of the anti-hermitian matrix X as: 

-X^ = X = -ZxX^E„ (3.1) 

where the matrix Ex is defined as: 

The Lie algebra of the matrices defined in Eq. (13.11 ) is isomorphic to the so{2M) algebra^. Since 
X belongs to a Lie algebra, it can be diagonalized: 

X = U"UU. (3.2) 

Here U is an element of a Lie group which is isomorphic to S0{2M). The diagonalization 
of Eq. (13.21) is equivalent to decomposing the X matrix in matrix polar coordinates. The radial 
part corresponds to the A matrix, A = diag (/A, —iX), where X is an M x M diagonal matrix of 
eigenvalues X = { Ai , . . . , Am}, and the angular part corresponds to the U matrix. 



B. Diagonal form 

Let us now consider the exponent of the Gaussian operator defined in Eq. (|2.7I) . After diago- 
nalization, we can write this operator expression in terms of the matrix X as: 

^zi(at.a)(u-Uu)(;). 

Since the canonical anticommutation relations are invariant under this transformation, we can 
define new fermionic operators 




(3.3) 



therefore, we can show that: 

G(H) =exp 



a 




and hence using anti-commutation relations it follows that 

G(A) = exp 

We can write Eq. (13.41) in normally ordered form using Eq. (12.141) . hence 



Ab^b-^Tr[A] 



G(A) = ^>-^Tr[A].exp 



(3.4) 



(3.5) 



C. Elementary properties 

We now wish to use this diagonal form to prove some elementary properties of these Gaussian 
operators. 

1. Positivity 



The Gaussian operator defined in Eq. (13.51) is defined in terms of real eigenvalues, Xj G 
(— oo^oo) ^ so that e'^j G (0,°°). Hence, if we consider the normally ordered form we notice that 



exp 



G{X) = e^2' 

= n^"''^'(l+b? (e^-l)^b,)>0. 



(3.6) 



Therefore the Gaussian operators with hermitian H matrices are themselves hermitian, positive 
definite operators in the fermionic Hilbert space. 



2. Normalization 



The Gaussian operators defined in Eq. (13.51) can be normalized to obtain operators A(H), such 
that Tr [A(H)] = 1. We first consider the trace of a single-mode case, which is given by: 



Tr 



: exp 



bUe^-l\b 



Tr 



On the other hand, including the exponential factor obtained during normal-ordering, 

Tr 



\+e^]e^^ =2coshf4 



(3.7) 



(3.8) 



Then for the single mode case we can introduce a normalized Gaussian operator A(A) such that 



Tr [A(A)] = 1, in the form A(A) = exp 



/[2cosh(A/2)]. 



For the general M-mode case we can transform back to the original fermionic operators, using 
invariance of the determinant under unitary transformations to obtain: 



A(H) 



e3 



iyl'Hy 



(3.9) 



det [2cosh(H/2)]' 

We note that this normalized form is identical to that used in Gaussian phase-space representations^"—. 



IV. GENERAL RESOLUTION OF UNITY 

A. General Gaussian Operators 

We wish to prove that the resolution of unity for the normalized hermitian Gaussian operator, 
A (H), defined in Eq. (13.91 ) is given by: 

J dllA{ll)P{lf) =1. (4.1) 

Here / is the fermionic identity operator, JH is the measure over hermitian matrices which 
conserves the nonstandard group symmetries of any of the four different symmetry classes de- 
fined by Altland and Zimbauer^. Here, ^ (H^) is a normalizable positive function of that 
is invariant under the transformation U, where U depends on the symmetry class. On taking a 
trace of the equation, we note from Eq (14.11 ). that the normalization of this function must satisfy: 

2-^fdnp{n^) = 1. 

Although our result is general, for definiteness we consider two options for the normalization 
function: 

/'(l)=dl)det[l+H2]"^ (4.2) 

p{2) ^ ^(2) ^jj2j j ^ (4_3) 



B. Matrix polar coordinates 

We now wish to use a unitary transformation to reduce the integral over matrix elements of H 
to matrix polar coordinates. The angular variables correspond to unitary transformations, while 
radial variables correspond to eigenvalues. The Jacobian for the transformation from the Cartesian 
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Table I. Indices of the different symmetry classes. 



Class 




a 


D 


2 





C 


2 


2 


DEI 


4 


1 


CI 


1 


1 



coordinates H = —iX to polar coordinates (A,U), is given by diagonalizing the matrix H defined 
in Eq. (|Z8l) , so that^ 



(4.4) 



Here dX = Y[f=\dXj, U is the transformation that diagonalizes the matrix H for the symmetry 
class under consideration, and A (A) is the Vandermonde determinant defined as: 



A(A)=A(Ai,...,Am)= n (A-Aj)=det 

l<i<j<M 



(4.5) 



The indices a and j8 of Eq. (|4.4I) . depend on the underlying nonstandard symmetry class. Their 
values are given in Table HI 

Eq. (14.11) in polar coordinates can now be written, after unitary transformation to diagonal 
operator form, as: 

jHA(H)p (h2) = I (vUv) I p (a') ^^ (a2) dX X 



M 

xn 



|Aj|"e : exp 



[2 cosh {Xj/2)] 



(4.6) 



Just as in the previous section, we have defined bj as a function of the unitary transformation U, 
according to Eq (13.31) . 



C. Evaluation of integrals 



We next wish to show that the result of the integral over radial variables A is a constant, inde- 
pendent of the transformation U. In order to evaluate the integral over the radial variables A , we 
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notice that the expansion of the operator term can be expressed as: 



Y[e : exp 



= ]^ +2sinh b]b 

Using this result to expand Eq. (14.61) we obtain: 

Mr M 

dXAp (a^) A(A)/' (a^) n 1^/ = / (^') (^') n 

;=i .7=1 

+ fdXp[x^)Ap (A2)n 

^ i= 1 



(4.7) 



[2cosh(A;/2)] 



I Ay I tanh 



A, 



.(4.8) 



Here we see that the integral over the operator terms b\ bi has tanh terms which are all odd in the 
Xj variable, while every other term is an even function of Aj. Therefore, from the parity of these 
functions, all the terms of the second integral in Eq. (14.81) vanish on integration over Aj. 

Since the operator terms are the only terms that depend on the unitary transformation U, it 
follows that the integral over the angular part is a constant, given by the relevant angular volume, 
which we will denote by C^. The value of depends on the corresponding non-standard sym- 
metry class under consideration, which defines the transformation U. The result of the integration 
over the eigenvalues then determines the normalization factor. 

Using the standard identity, exp (Aj/2) = cosh (Aj/2) — sinh {Xj/2), Eq. (14.81) can be written 

as: 



JHA(H)P(h2) =C^2-^ I dXp{X?') aP (a') n|A, 

=2-^cu r dXp{x^)A^ (^')ni^,/r- 



A/ 

1 - tanh ( 



(4.9) 



In order to obtain Eq. (14.91) . we have used the result that the integral over the tanh (Aj/2) terms 
vanishes, because clearly tanh (Aj/2) is an odd function of Aj while the other terms are even. 
Next, we recall that the original definition of P was such that it was normalized. Therefore, using 
matrix polar coordinates to evaluate Eq (|4.1I) . we find that: 



JHP (H^) 



dXP^k'^A^ (^'jni^i-r- 



(4.10) 
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For any non-standard symmetry class it is therefore possible to express the resolution of unity 

as: 



JHA(H)/'(h2) 

dXp(^X^)AP^?i^)fim". (4.11) 

The value of both the angular volume and the radial integral will depend on the corresponding 
non-standard symmetry class under consideration. For the radial part, the integral will also depend 
on the different choices of the normalization function P{X ). We notice that as long as we can 
perform the integration over the angular and radial part, then it is possible to obtain a resolution 
of unity for any of the non-standard symmetry classes. An explicit result for these integrals in the 
case of class D symmetry is given next. 

We note there is an important consequence of this result in random matrix theory. Suppose we 
consider a statistical random matrix mixture of finite temperature canonical ensembles given by: 

p = J dUP (H^) exp [-/3H] oc /, (4.12) 

where H is the linearized Bogoliubov-de Genne Hamiltonian given by Eq (12.91) , and JH is a mea- 
sure over one of the four nonstandard symmetry classes. It follows that any correlation function 
or moment evaluated in this ensemble is simply an average over the identity operator, independent 
of temperature, symmetry class or the details of the distribution P (H^) . 



D. Symmetry Class D 

In order to give the values of the angular volume and radial integrals defined in Eq. (14.1 II) . 
we will consider the most general symmetry class D. That is, we are going to integrate over the 
transformations described in Sec. IIII A[ We note, however, that similar results in any of the non- 
standard symmetry classes can be found. The main differences are in the values of the integration 
volumes, and in convergence properties which depend on the details of the Jacobian and matrix 
polar coordinates. 

For symmetry class D, which is the largest symmetry class, the values of the indices a and /3, 
given in Table H are and 2 respectively^. In this case the value of the integral over the angular 
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part, CU = VUV is given below. Details are given in the Appendix: 

Hence, for the non-standard symmetry class D, the integral defined in Eq. (14.111) is: 

1 = J JHA(H)/'(h2) 

/oo 
/'(a2)a2(^A')jA. (4.14) 

1. Determinant normalization 



In order to evaluate the integral over the radial variable X , we have to consider one of the two 
options for the normalization P (^X,^^. We first consider the normalization by the determinant 
given in Eq. (|4.2I) . Hence the integral over the variables X is: 



/oo poo lYi 

jaa2 (a^) (a^) = c« / dk^^{x^)Y[{l+ljy^'. (4.15) 
-OO iV — oo J — 

The integral of the right hand side of Eq. (|4.15l) is known as Selberg's Integral-^^. Selberg's 
formula is given below. It is valid for integer n and complex a, /3, y with Rett > 0, Re/3 > and 
Rer>-min(i,Sf,g^): 

poo poo /7 

Jo Jo ,=1 



w-1 



r(l + y+jy)r(a + jy)r(/3 + jy) 
U rfl+y)r(a + /3 + (n + j-l)y) ■ ^ 



1 



Defining Vj = Xj, dv = iXdXj, dXj = jV ^dv, and setting n = M, a = 1/2, y = 1 and -a - 
/3 — 2y(n — 1) = —2p, so that /3 = — 2M-|-2j9-f 3/2, we require p> M which gives the result: 

Therefore, one fermionic resolution of unity in the most general symmetry class D is given by: 

J JHA(H)/'(i) (H^) , (4.18) 

where the normalization constant C^^^ is given by: 



^M(M-i) 11 r(-2M + 2;? + j + |)' 
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2. Gaussian normalization 



Now we consider the second option for the normalization function, given in Eq. (14.31) . where 
we normalize the distribution by a Gaussian c-number function of the eigenvalues, as often used 
in random matrix theory. Since this normalization is another even function, the integral over 
Gaussian operators reduces to a term proportional to the identity operator, as before. We now wish 
to evaluate the normalization constant, to obtain a resolution of unity. In this case the integral over 
the radial part is: 

/oo poo 
dXA^ ^(2) (a^) = C(2) / dXA^ (a^) e-^'^^^^J. (4.20) 

The integral of Eq. (14.201) is a Selberg type integral related to the Laguerre polynomials^: 

I A (^)rtl exp (4) * = 2«-M«-.) n r(i+;r)r(«+rO-i)), 

J-oo J-oo ,=i \ ^ / ;=i r(l + 7) 



We set a = 1/2, 7 = l,n = M and JCy = 4-pXj. Hence we obtain: 

jy^) (a^) a2 (a^) ja = (2;.)-^^-^) nr(i + j)r (j - [ 

In this case the resolution of unity is 

/ = J JHA(H)/'(2) (h2) , (4.21) 
where the normalizing constant is given by 

\M(M-l/2) 

d2)=2^M^j . (4.22) 

V. NUMBER-CONSERVING GAUSSIAN OPERATORS 

By analogy to the general Gaussian operators, we wish to investigate if there is a similar ex- 
pression for the fermionic identity operator in terms of the number-conserving Gaussian operators. 
These operators, Ga^ (h) , defined in Eq. (12.151) can be normalized so that Tr [GAr(h)] = 1. The 
normalized number-conserving Gaussian operator An (h) is: 



iTr(h) 



^^^^^^ = det(2cosh(h/2)) ■ 
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a' 



gh-i 



a 



(5.1) 



In section HTCl we introduced the variable u = e^. We notice that if we define the variable u as 



u = n 



I, we obtain the expression for the normalized number-conserving Gaussian operator 



defined in terms of the stochastic Green's functions^ii^ for particles, n and holes ii = I — n: 



Ayv(h) = det [n] : exp a"^ [n.-^ -21] a 



(5.2) 



A. Matrix polar coordinates 

We wish to investigate if there is a normalization factor ^(h) that generates a resolution of 
unity for the normalized number-conserving Gaussian operators, Aa? (h) defined in Eq. (15.11) . of 
the form: 

J AN{h)^{h)dh = I. (5.3) 

Here dh is an integration measure over the hermitian matrices h and ^(h) is a normalization 
function. This normalization function is defined in order to ensure the convergence of the integral. 
We can consider, for example, the following option for the normalization, in analogy with the 
previous case: 

C 



(5.4) 



det [cosh (h/2)]^' 
where C is a normalization constant and we require p > I. 

Hermitian matrices can be decomposed in polar coordinates as^ h = UAU^, where U is an 
unitary matrix and corresponds to the angular coordinates, while A corresponds to the radial coor- 
dinates and is a diagonal matrix, X = {Ai, . . . , Am}. The Jacobian of the transformation from the 
cartesian coordinates h to polar coordinates (A, U) is given by^: 



dh = dVA^{X)Y[dXi. 



(5.5) 



Here dV is the normalized Haar measure over unitary matrices and A(A) is the Vandermonde 
determinant defined in Eq. (14.51) . 

B. Evaluation of integrals 



Hence, on diagonalizing the Gaussian operator: 
An (h) = : exp 



'1^2 cosh (A;/2) 



(5.6) 
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where we have defined b = V^a. Therefore Eq. (|5.3I) in polar coordinates is: 



Aw(h)^(h)Jh = C / dVA^{X)Yl 

7=1 



M 



: exp 



2cosh(A;/2) 



d?ij. 



We will now focus on the integral over the radial part. Using the expression of the operator term 
defined in Eq. (14.71) . we obtain: 



-Xj/2 



An (h) ^{h)dh = C J dUA^ (A) ]^ (^1 +b] (^e^j - l) bj 

= c/jUA2(A)n TT+ +1 



2cosh(Aj/2)^+^ 
sinh (A;/2) 



;=i \cosh (Aj/2)'' cosh(A//2)' 
Here we notice that the Vandermonde determinant term is not an even function of each eigen- 
value. Hence, we do not obtain a resolution of unity following simple parity arguments with an 
even weight function, as in the previous section. The reason is that when expanding the terms of 
the Vandermonde determinant, !!;< j {^i ~ ^j) i for different values of i and j, we eventually obtain 

/V 4. /V 

operator terms bjbj that depend on even and odd functions of Aj for a fixed value of p. 

On the other hand, if ^(h) is not unitarily invariant, and includes terms that when multiplied 
by the Vandermonde determinant give an even function of A/, then we reach a different conclusion. 
An example of this is if the weight function has the form: 



^(h) = 

In this case, it is clear that: 

A2(A)^(h) 



n i^i+^jf 

\<i<i<M 



n (^-^ y 

l<i<j<M 



(5.7) 



(5.8) 



With such a weight, the integral corresponding to the operator terms in Eq. (15.31) vanishes, 
because of the parity of the functions: the function sinh (Aj/2) is odd while the other are even 
functions of Ay. In this case we can obtain a resolution of unity just as in the nonstandard symmetry 
case. Therefore, for the number-conserving Gaussian operators, we can obtain a resolution of 
unity if the normalization factor cancels the parity violation of the Vandermonde terms. Thus, we 
have shown that an expansion of the identity operator is also possible with the number-conserving 
Gaussians. As this is not a simple trace or determinant, we cannot easily express this in terms of 
the original coordinates h. These arguments obviously do not exclude other routes to obtaining a 
resolution of unity with these operators. 
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VI. SUMMARY 



We have studied a subset of the general fermionic Gaussian operators, the hermitian positive 
definite fermionic Gaussian operators, which belong to a non-standard symmetry class of random 
matrices. We have proved that there are simple resolutions of unity for these operators, for any of 
the non-standard symmetry classes, and for any integrable even distribution of eigenvalues. This 
resolution of unity is defined for the entire Hilbert space, not for a coset or subset space as in the 
case of the fermionic coherent states. 

Our proof is based on considering the hermitian Gaussian operators as an operator basis de- 
pending on a continuous hermitian matrix. We use polar coordinates for skew- symmetric matri- 
ces, which lead to an expression for the resolution of unity in terms of integrals over eigenvalues. 
Our result appears similar to related expressions in random matrix theory. In order to obtain this 
unitarily-invariant expression for the resolution of the unity, the normalizing factor must be an even 
function of the eigenvalues. This suggests that distributions of random matrices in the nonstandard 
symmetry classes, if used to average over canonical density matrices, should generally have other 
types of distribution. Otherwise the physical behaviour will simply correspond to a unit density 
matrix. While this is not impossible, it is unlikely to be generic. In the case of the symmetry class 
D, we give the values of the constants corresponding to the angular and radial integrals, with two 
different options for the normalizing factor. 

These resolutions of unity for the Gaussian operators can be used to derive mathematical iden- 
tities for other physical applications. As an example, our results can be used to construct a positive 
fermionic distribution or fermionic Q-function. This will be carried out elsewhere. 
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Appendix A: integral over the angular coordinates 

In this Appendix we evaluate the integral over the angular coordinates U, whose expression is 
given in Eq. (|4.13l) . for the non-standard symmetry class D. 
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We wish to evaluate the constant = J (U^JU) . This is the integral over the angular variables 
U. It is the angular part of the Jacobian of the matrix transformation from Cartesian coordinates 
H = iX to polar coordinates (U, A): 



J dHP (h2) = J (U^JU) J dXA" (a^) flP{X}) . (Al) 

The value of the integral over the angular variables U is given by the ratio of an integral over 
cartesian coordinates H, to an integral over radial coordinates A : 

Since the constant is evaluated as a ratio, we will perform the calculations using the Gaussian 
form of the normalization function. This gives simple integrals in Cartesian matrix coordinates. 
For these calculations, we therefore use: 

P (h2) = exp [-pTr [H^] ] . (A3) 

Consequently, we now have to evaluate: 

CU ^ /^Hexp[-pTr[H^]] ^^^^ 

Radial integrals 

The radial integral over the eigenvalues Xj in Eq. (|A4I) is a Selberg type integral related to the 
Laguerre polynomials^: 

In our case, we consider 7 = 1, a = 5, so that the radial integration gives: 

fjXA^ (A^) e-'^^^^'f = (2;.)-^(^-i) nr(l + j)r (j - ^) . (A6) 

Cartesian integrals 

The next step is evaluate the integral over the 2M x 2M hermitian matrix H, defined as: 

( h A \ 

H= 1 , (A7) 
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with hap = h*p^ and Ao,p = —Ap^,. We consider the integral over the matrices h and A so that: 

r M 

JHexp[-rr(pH2)] = / ]^ J/z^^J/z^^.JAf/Aj^-exp [-pTr (H^)^ . (A8) 

•J i=l i<j 



The trace of can be written as: 

'J 



M M 

E 



=2£/.^+4j: {h^j)'+{h^j) 

i 

Hence, we can use the result that: 



. 2 

V 



to obtain the overall Cartesian integral of: 

-N M(2M-\) 



where we have integrated over both the diagonal and the off-diagonal terms. 



Normalization constant 



Using the results of Eq. (IA6I) and Eq. (lAllI) we obtain: 

M-l 



(A9) 



dhe-^P^'' = J^, (AlO) 
V 



j JHexp [-pTr{n^)] = ^ ^-MiM-i) ^ ^^^^^ 



|ir(2+j)r(j + i) 
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